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The two inequivalent valleys in graphene preclude the protection against inter-valley scattering 
offered by an odd-number of Dirac cones characteristic of Z 2 topological insulator phases. Here 
we propose a way to engineer a chiral single-valley metallic phase with quadratic crossover in a 
honeycomb lattice through tailored VSN x or 3N x 3N super lattices. The possibility of 

tuning valley-polarization via pseudo-Zeeman field and the emergence of Dresselhaus-type valley- 
orbit coupling are proposed in adatom decorated graphene super lattices. Such valley manipulation 
mechanisms and metallic phase can also find applications in honeycomb photonic crystals. 

PACS numbers: 68.65.Cd 71.10.Pm 73.22.Pr 73.43.Cd, 


Introduction —. Honeycomb Dirac materials have two¬ 
fold de^nerate band structures with inequivalent KK' 
valleys [iMsl, whose origin can be traced back to the bi¬ 
partite nature of honeycomb lattices (A and B triangular 
sublattices). This binary valley degree of freedom has 
led to proposals of valleytronics applications iini that 
leverage the valley pseudospins in a manner analogous 
to electron spins in spintronics applications. A distinct 
scenario is that of single (odd-number) Dirac-cone in Z 2 
topological insulators 0 where their surface states are 
effectively decoupled from each other due to their distant 
spatial separation. Therefore, a single Dirac-cone struc¬ 
ture is desirable when we require a Hamiltonian that em¬ 
bodies the chiral anomaly of Dirac fermions [lH and at 
the same time is protected against inter-valley scattering. 

In this Letter, we propose to engineer a single valley 
phase in 2D honeycomb Dirac materials through \/3N x 
V3N or 3N X 3N superlattices that fold and couple the 
inequivalent KK' valleys into the same T point. We show 
that the corresponding effective Hamiltonians for top-site 
adsorbed superlattices exhibit uniform inter-valley cou¬ 
pling and valley-orbit coupling mechanisms that resemble 
the conventional in-plane Zeeman fields and Dresselhaus 
spin-orbit coupling of the electron spins [ESliilii. The 
pseudo-Zeeman field and pseudospin-orbit coupling allow 
to control valley polarization coherently, while the latter 
one further indicates the possibility of controlling valley 
polarization via electric fields. Moreover, together with 
the coexisting sublattice potentials, we find that inter¬ 
valley coupling can drive a topological phase transition 
from a quantum valley-Hall phase into a chiral single¬ 



FIG. 1: (color online) Schematic representation of inter¬ 

valley coupling adatom superlattices and their respective Bril- 
louin zones, (a) and (b) are respectively primitive and recip¬ 
rocal lattices for the top adsorption in \/3 x \/3 graphene su¬ 
percells. The red lines represent the Brillouin zone of pristine 
graphene. 


valley metallic phase with quadratic band crossover. We 
also propose that such inter-valley coupling mechanism 
and metallic phase can be observed in photonic crystals. 

Inter-Valley Coupling —. When the \/3N x V3N or 
3N X 3N supercells are tailored on a honeycomb lattice, 
KK' valleys couple and fold into the T point, as illus¬ 
trated schematically in Fig. [Hb) showing the reciprocal 
lattices for both 1x1 (red) and x ^/3 (black) super¬ 
cells. For definiteness, here we only focus on the top-site 
adsorption as shown in Fig.^a) and leave the discussion 
of the effective Hamiltonians for bridge- and hollow-site 
adsorption in the Supplemental Material. For top-site 
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adsorption in a \/3 x ^/3 supercell, the six atoms in each 
primitive cell can be classified into three different cate¬ 
gories: (i) one at the adatom site, (ii) three at the nearest 
neighbor sites, and (iii) two at next-nearest sites. We rep¬ 
resent the corresponding site energies as ui, U2, and 1^3, 
and set 1^3 = 0 as the reference value. Assuming that the 
adsorption sites belong to sublattice “A”, the real-space 
tight-binding Hamiltonian can be written as: 

Ht = Ho + api + ^ bl_^_gbi+s, (1) 

i i S 

where runs over all adatom sites. Here Hq = 
+ h.c.) is the band Hamiltonian with to 

being the nearest-neighbor hopping energy, and a| (6|) 
is the creation operator of an electron at i-th A(B) site. 

The Brillouin zone of pristine graphene can be rep¬ 
resented through three copies of ^/3 x \/3 graphene 
supercell’s Brillouin zone as displayed in Fig. mb). 
By denoting j-th {j=l- 3 ) center as Kj, the operator 
ai can be expanded in momentum space as: = 

7^ J 2 j exp[-i(Kj + k) ■ Ri]aj,jk, where Nq is a nor- 
malization factor, and k runs over the Brillouin zone of 
X graphene supercell. The Hamiltonian of Eq. ([6]) 
in momentum space is: 

H,ik) = Hoik) + ( 2 ) 

jj' 

where Ho{k) = -toJ 2 j{Xjkaj,.bj^k + h,.c.) describes 
the kinetic energy of pristine graphene with Xjk = 
and . The Kj 

(j=l- 3 ) are respectively wavevectors of K, K', and T 
points. The last two terms give sublattice potentials 
when j = j' which are different for AB sublattices due to 
inversion symmetry breaking. When j ^ j', they give rise 
to inter-valley coupling through a finite ui contribution 
while U2 contribution vanishes due to the phase interfer¬ 
ence (^KK'= 0 ). By block diagonalization, the low-energy 
effective Hamiltonian can be further obtained: 

= Uo^VF{kx(Jx^Tzky(jy) + ( 3 ) 
Ao 

+ 

where Uo={A2 + U2)/2 and Ai = (A2 — U2)/2 with A2 = 
1^1/ 3 . The third term reflects the effective potential im¬ 
balance through a mass term of magnitude Ai and the 
last term describes inter-valley coupling through the Tx 
operator. We note that the coupling between K and K' 
valleys only occurs at “A” sublattice with the coupling 
amplitude A2 depending on ui linearly. Such an inter¬ 
valley coupling acts on the valley pseudospin as an effec¬ 
tive Zeeman field that can be used to control the valley 
polarization coherently in valleytronics devices. 

When the nearest neighbor hopping terms of superlat¬ 
tice Hamiltonians are allowed to change hy St = t — to 


due to the influence of the adatoms, the real-space tight- 
binding Hamiltonian in Eq. ^ acquires an additional 
term H' = Yl\i where the index i runs 

over “A” sites right underneath the adatoms and the j 
sites represent the three nearest “B” sites [2^ . The mod¬ 
ified effective Hamiltonian becomes: 

Hf^'ik) = U'o+v'pia^lrkx + (JyT^ky) + 

A' 

+ ^(1<T + Cr^)ra; V5(Ta;(ra;fca; - (4) 

where ([/q, A'l) have same forms as (t/o, Ai) by changing 
A2 to be A2 = 3 i 4 ito/(t-l- 2 to)^, and the Eermi velocity is 
modified to be v'p = VF{ 2 t to)/{t 2 to ) • The last term 
in Eq. o can be identified as a Dresselhaus-type valley- 
orbit interaction of strength vs = VF{t — to)/{t ^ 2 to) 
coupled with a sublattice-fiip potential. This term also 
couples different valleys and implies the possibility of ma¬ 
nipulating the valley degree of freedom by external elec¬ 
tric field in a manner analogous to the control of electron 
spin by electrical means via spin-orbit coupling. 

Single- Valley Metallie Phase —. Adatom superlattices 
lead to both inter-valley coupling and inversion symme¬ 
try breaking potentials, and it is easy to understand 
that each term can independently contribute in open¬ 
ing a Dirac point gap when they are viewed as uniform 
in-plane xy and 2: contributions to the pseudospin fields 
in the Dirac Hamiltonian , where the former shifts 
the position of the Dirac points in momentum space and 
the latter introduces an inversion symmetry breaking gap 
in the Dirac cone. Here we show that when those ef¬ 
fects are present in a superlattice, a topologically dis¬ 
tinct single-valley phase can be engineered. We begin 
considering for sake of clarity the top-adsorption config¬ 
urations neglecting the modification of the hopping en¬ 
ergy in the band Hamiltonian and setting the site en¬ 
ergies at all “B” sublattices to assume a constant value 
(z.e., Ub = U2 < 0 ). When = 0 , the site energies at all 
“A” sublattices are identical, z.e., Ua = 0 . This leads to 
vanishing inter-valley scattering and the imbalanced sub¬ 
lattice potentials open a quantum valley-Hall gap at the 
Dirac points, where the doubly-degenerate massive Dirac 
cones are folded as a single valley around the T point but 
remain distinguishable [See Eig.[ 2 ][a)]. When we allow ui 
to take negative values, we find a gradual decrease of the 
inversion symmetry breaking induced gap | Ai | and an in¬ 
crease of inter-valley coupling strength IA2I that lifts the 
degeneracy of the conduction bands splitting by a magni¬ 
tude of 2A2 [See Eig.[ 2 ][b)]. The simultaneous presence of 
both terms breaks the particle-hole symmetry and leads 
to a smaller bulk gap A' = |2Ai + A2I. 

When ui is even further decreased and reaches a crit¬ 
ical value of ui = 3i42/2, the bulk gap A' completely 
closes. As shown in Eig. [2][c), we achieve a single band 
touching point at T formed by a Dirac-cone centered at 
the edge of the parabolic valence band. In this limit 
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(a) (b) (c) (d) 



FIG. 2: (color online) Topological transition from a quantum 
valley-Hall insulator to a single-valley phase as a function of 
the parameters ui. Here, we set U 2 to be fixed with 1^2 < 0 
and 2/3 = 0 . — 2 A 2 corresponds to the local band gap from 
the inter-valley scattering. A' measures the bulk (local) band 
gap from the competition between inter-valley coupling and 
sublattice potentials. The progressive decrease of 2/1 leads to 
a complete closure of the quantum valley-Hall gap and then 
transitions to the single valley phase by reversing A'. 


where the bulk gap is closed, the valley-Hall effect is ab¬ 
sent and the valleys are no longer distinguishable. When 
we allow even smaller values of 2 / 1 , the inter-valley cou¬ 
pling strength IA 2 I further increases, while the magni¬ 
tude of the staggered sublattice potentials |Ai| first de¬ 
creases to zero then increases again [ 2 ^ . When the inter¬ 
valley coupling is strong enough, a valley-mixed metal¬ 
lic phase with quadratic band crossover is engineered as 
displayed in Fig. [Sfd). In this limit the edges of the 
lower energy bands are distant from the corssing point 


by A' 


2|Ai| - IA 2 I 


[see Supplemental Material for 


details of the tight-binding band structure]. 

A detaild analysis of the low energy bands reveals that 
the parabolic dispersion at F point is chiral and formally 
similar to the band dispersion near K/K' valley in Bernal- 
stacking bilayer graphene. When IA 2 I ^ |^i|: the low- 
energy Hamiltonian of the quadratic touching bands can 
be further simplified as: 


Hf (fe) 


U^' + ak^ 




0 

TT^ 0 


(5) 


which is represented on the basis of “B” sublattice 
from both K and K' valleys. Here, we define = 
Uo - Ai, e = kl^ kl, a = Ai4/(A2 - A?), and 
P = A 22 ;|./(A 2 — A^). The last term couples states be¬ 
tween valleys K and K' with n = kx ^ iky ^ and gives rise 
to the quadratically dispersing Fermi point band struc¬ 
ture. Such a two-fold degeneracy at the crossing point, 
which also appears in Figs. [2][a)-|2][c), is protected by the 
Csv symmetry, since these two basis functions form a 
two-dimensional irreducible representation of the corre¬ 
sponding point group. For other adsorption geometries, 
e.g., hollow- and bridge-adsorption with respectively C^y 


TABLE I: Inter-valley coupling mechanisms for different ad¬ 
sorption geometries. 


Adsorption Site 

Symmetry 

Inter-valley Coupling 

Top 

Csv 

(1 + az)rx 

Hollow 

Cqv 

rx^y 

Bridge 

C2V 

r’cc 1(T 


and C 2 V symmetries, the band structures as well as inter¬ 
valley coupling mechanisms become completely different 
as listed in Table 

The main difference of this band crossover from the 
case of bilayer graphene is that here we have only a single 
Dirac parabolic dispersion. This is of interest, because 
it provides an ideal platform to study the single Dirac- 
cone transport phenomena of Z 2 topological insulators 
and allows to explore the chiral anomaly of single val¬ 
ley physics that is not compensated by its time-reversal 
counterpart. For example, if broken symmetry gapped 
phases are developed in the presence of electron-electron 
interactions l22l424| . a mass sign dependent spontaneous 


orbital moments will develop per spin-valley 23-2S|. In 
our single valley phase, it is expected that when the Fermi 
surface lies at the crossing point, a quantum anomalous 
Hall ground state will develop when both spin compo¬ 
nents have the same mass, or alternatively a quantum 
spin-Hall state will be present when the masses for each 
spin term have opposite signs [l^. Besides, a supercon¬ 
ducting phase can also be expected when the Fermi sur¬ 
face is shifted away from the crossing point [HI . Whereas 
the energetically favored ground state depends on de¬ 
tails of the band Hamiltonian and the models for the 
electron-electron interaction, further control of quantum 
phase transitions should be achievable by means of exter¬ 
nal magnetic fields coupling with the spontaneous orbital 
moments. Furthermore, in bilayer graphene, the magni¬ 
tude of the gaps predicted in a Hartree-Fock theory with¬ 
out d yna mical screening is on the order of a few tens of 
meV [ 2 ^ whereas experimental gaps turned out to be an 
order of magnitude smaller ^2 meV [ 2 ^ due to the ex¬ 
ponentially increasing screening feedback when the gaps 
are small. Thus, it is expected that substantially larger 
gaps can develop, if flatter bands can be tailored when 
the leading parabolic dispersion coefficients can be made 
smaller than the one used in bilayer graphene. Moreover, 
in the presence of strong magnetic field, the anomalous 
Landau-level quantization can also be expected as that 
in bilayer graphene case [30| . 

Photonic-crystal bands —. Experimental realizations 
of periodic graphene superlattices could take advantage 
of substrates that can generate the 3 x 3- or V3 X ^/3- 
type superstructure, like EuO(lll) M and Ag(lll) sub¬ 
strates [111. There are also other methods for engi¬ 
neering such kind of superstructures, e.^., silicene on 
Ag(lll) substrate [HI, InSb(lll) surface Q, artificial 
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FIG. 3: (color online) Upper panel: Schematic representa¬ 
tion of honeycomb photonic crystals with a VS x Vs peri¬ 
odicity. di and d 2 denote the primitive vectors. The dis¬ 
tance between nearest columns is set to be a and the slab 
width is d = 0.1a. (z = 1-3) label the radii of differ¬ 
ent columns. In our simulation, each column is chosen to 
be inhnitely long. Lower panel: Photonic band structures of 
transverse-magnetic modes along high-symmetry lines for dif¬ 
ferent radii ri = 0.18a (a), 0.23a (b), 0.28a (c), and 0.32a (d), 
respectively. Here, we set r 2 = 0.25a and rs = 0.18a. 


organic molecular lattice [35|, or patterned two dimen¬ 
sional electron gas with well-established experimental 
technique [s^. The applicability of our theory depends 
on the degree of the achievable commensurability with 
the crystal structure of honeycomb lattices. It is note¬ 
worthy that, since our model is spin independent, it can 
also apply to Bosonic systems like cold atoms or pho¬ 
tonic crystals Q in honeycomb super lattices. One pos¬ 
sibility is to use honeycomb photonic crystals made of 
silicon columns linked by thin silicon slabs as shown in 
the upper panel of Fig.[3l and use electromagnetic waves 
with transverse-magnetic modes in the xy plane. The 
corresponding site potentials and hopping energies for 
the photonic crystal setup can be controlled through the 
column radius r and the link width d. The confinement 
radii allow to tune the concentration of electrical-field 
energy of the harmonic modes. 

If the columns’ radii are identical and the connecting 
slabs have the same width, the two dimensional pho¬ 
tonic band structure for transverse-magnetic modes [1^ 
obtained from the finite elements method Hi shows 
two linearly dispersing Dirac cones, clos^ resembling 
the band structure of pristine graphene 0 , To 

model the a/3 x graphene supercell, we first classify 
the columns’ radii into three categories Vi (i=l-3) with 
rs = 0.18a as a reference, and the link width is chosen to 
be d = 0.10a with a being the distance between two near¬ 
est columns. Figures [3l^a)-|3l^d) display the photonic band 


structures for different ri = 0.18a (a), 0.23a (b), 0.28a 
(c), and 0.32a (d) at fixed r 2 = 0.25a along high symme¬ 
try lines. One can observe a topological phase transition 
from an insulator to a single-valley metallic phase when 
ri is progressively increased [See the highlighted regions] 
in a way closely similar to the behavior of the electronic 
band structure shown in Fig. O 

Discussions and Conclusions —. We presented the the¬ 
ory for the inter-valley coupling mechanisms in \/3 x \/3 
graphene super cells that act as in-plane pseudo-Zeeman 
fields or pseudospin-orbit coupling. Both contributions 
can be used to tailor valley pseudospins of honeycomb 
lattices. Especially, the Dresselhauls-type valley-orbit 
coupling makes it possible to control the valley polar¬ 
ization via electric means. These valley coupling mecha¬ 
nisms have important implications in valleytronics, where 
the coherent control of valley polarization is yet a grand 
challenge due to the missing counterpart mechanisms 
of spin-orbit couplings or magnetic fields for spintron- 
ics. Moreover, our theory also suggests strategies for 
engineering single-valley electronic structure in conven¬ 
tional Dirac materials with two inequivalent degener¬ 
ate valleys by folding them together. The single-valley 
phase can be manipulated by combining inter-valley cou¬ 
plings and imbalanced sublattice potentials originated 
from the inversion-symmetry breaking. By increasing the 
strength of inter-valley coupling from zero, a topological 
phase transition can take place from the quantum valley- 
Hall phase to a chiral single-valley metallic phase with 
quadratic band crossover that resemble the electronic 
structure of a half Bernal-stacked bilayer graphene. A 
concrete proposal for such a single-valley phase is pre¬ 
sented in honeycomb photonic crystals. We verified that 
all these findings can also be realized in 3 x 3 honeycomb 
supercells. 
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SUPPLEMENTAL MATERIALS 
INTRODUCTION 

In the main text, we presented the theoretical formula¬ 
tion for engineering a single-valley phase by folding K and 
K' valleys of graphene into the F point via the top-site 
adsorption in a \/3 x \/3 supercell. In the following we 
supplement the information in the main text by providing 
more explicit details to facilitate understanding. In Sec¬ 
tion II, we provide a detailed analysis of the phase tran¬ 
sition by exactly diagonalizing the tight-binding Hamil¬ 
tonian, where an interesting partner switching of the en¬ 
ergy bands occurs along with the phase transition from 
quantum valley Hall insulator to a metallic phase wit 
quadratic band crossover. In Section III, we present the 
details about the calculation of the photonic band struc¬ 
ture and show the Dirac cones of honeycomb photonic 
crystals with uniform radii of slabs. In the main text, 
we have only focused on the top adsorption that breaks 
the inversion symmetry. When the atoms are adsorbed 
at hollow or bridge sites, the inversion symmetry is pre¬ 
served. For completeness, in Section IV, we present the 
physical origins of inter-valley couplings for hollow- and 
bridge-site adsorption and the corresponding band struc¬ 
tures. 


TOP ADSORPTION IN ^3 x ^3 GRAPHENE 
SUPERCELL 

For the top adsorption, we choose x supercell 
as an example [See Fig. 311^)]. In this superlattice, if 
we only consider the site energy modification from the 
adatoms, all these site potentials Ui can be classified 
into three types: ui for the sites right underneath the 
adatoms, U 2 for the sites nearest to adatoms, and for 
all other sites. For simplicity, 1^3 = 0 is chosen as a refer¬ 
ence. Therefore, the 7 r-orbital tight-binding Hamiltonian 
can be written as: 

iJt = -ffo + WlX “b* + “ 2 X X (6) 

i id 
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FIG. 4: (color online) (a): Primitive cell for the \/3 x \/3 
graphene supercell with top-site adsorption, m {i = 1-3) rep¬ 
resent the site energies of the three type atoms, (b): Recipro¬ 
cal lattices for 1 x 1 (in red) and x \/3 (in black) graphene 
supercells. 



FIG. 5: (color online) Evolution of tight-binding band struc¬ 
tures for the \/3 x VS graphene supercell for different ui = 0 
(a), 81 ^ 2/4 (b), 81 ^ 2/2 (c), 7 U 2 I 2 (d) at fixed U 2 < 0. Note: 
Bands’ partner-switching happens in the last panel after the 
topological phase transition. 


where runs over all adatom sites. Here Hq = 
+ h.c.) is the band Hamiltonian with to 

being the nearest-neighbor hopping energy, and a] {b]) 
is the creation operator of an electron at i-th A(B) site. 
The second term is the site energy induced by adatoms 
possessing the translation symmetry of VS x ^/3 super¬ 
cell. 

In the corresponding momentum space, the Brillouin 
zone of pristine graphene (red lines) is three times larger 
than that of the ^/3 x \/3 supercell (black lines) are plot¬ 
ted in Fig. ^ll^b). Thus, one can divide the Brillouin zone 
of pristine graphene into three parts with centers Ki, 
K 2 , and K 3 being wavevectors of K, K', and F points 
of graphene, respectively. Each part is the same as the 
Brillouin zone of x graphene supercell. Therefore, 
the operator can be expanded as: 

a* = exp[-i(Kj + k) ■ Ri]aj,fe, (7) 

h j 

where k runs over the Brillouin zone of \/3 x ^/3 supercell. 
By substituting Ect. (j7J) into Ect. m, the tight-binding 



EIG. 6 : (color online) Schematic of the staggered sublattice 
potential 2 |Ai| and inter-valley scattering IA 2 I as function of 

Ul. 


Hamiltonian in momentum space can be expressed as fol¬ 
lows: 

Rt(fe) = Hoik) + (8) 

jj' 

where Ho{k) = -to ^j{XjkCi]^k^j,k + h.c.) is the kinetic 
energy of pristine graphene with Xjk = , 

and ^jj' = xhe last two terms give sub¬ 

lattice potentials when j = j' which are different for AB 
sublattices due to inversion symmetry breaking. When 
j ^ /, they furnish inter-valley coupling where ui con¬ 
tribution is finite while U 2 contribution vanishes due to 
the phase interference (^kk'= 0 ). 

By exactly diagonalizing Eq. m, the tight-binding 
band structures for \/3 x ^/3 supercells with a negative U 2 
is plotted as shown in Fig. different ui. One finds 

that when ui decreases from zero, the band gap from im¬ 
balanced sublattice potentials [See Fig. 3Sl[a)] gradually 
decreases due to the emergence of inter-valley coupling 
[See Fig. 3Hb)]. When ui further decreases to a crit¬ 
ical point, the bulk band gap completely closes and a 
linearly dispersed Dirac-cone is formed centered at the 
edge of the parabolic dispersed valence band. For even 
smaller however, no bulk band gap reopens while a 
local band gap opens due to the increasing of inter-valley 
coupling. Such a local gap lifts the three-fold degenerate 
at F point and form a chiral single-valley metallic phase 
with quadratic band crossover. Interestingly, along with 
the phase transition, the third and fourth energy bands 
represented by green and purple lines switch partners. 

To understand the phase transition process more 
clearly, we first obtain the effective Hamiltonian by uti¬ 
lizing Lodwin block diagonalization method. In this 
method, we first rearrange the Hamiltonian in Eq. m to 
divide it into two blocks: one 4x4 matrix representing 
low energy bands contributed from valleys KK’, and the 
other 2x2 matrix representing high energy bands from F 
valley of graphene. Their coupling is much smaller than 
the hopping energy to and thus, the influence of high 
energy bands on the low energy ones can be obtained 
by perturbation method. Through a block diagonaliza- 
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FIG. 7: (color online) Left panel: Energy levels at F point 
for ui = 0. Both energy levels are double degenerate. Right 
panel: Energy levels at F point for ui = 81 ^ 2 /4. The degen¬ 
erate of the upper energy level in left panel is lifted while the 
lower one remains degenerate. 

tion 0, the effective Hamiltonian can be obtained as 
below: 

= Co + + AiCr^ (9) 

Ao 

where Uq = (A 2 +r^ 2)/2 with A 2 = and Ai = (A 2 — 
u^j2. The third term reflects the imbalanced sublattice 
potentials, while the last term describes the inter-valley 
coupling occurring at sublattice “A”. 

To show the dependences of inter-valley coupling and 
imbalanced sublattice potentials on the magnitudes 
of 2Ai and A 2 are plotted in Fig. 3S1 It shows that 
both terms depend on linearly. When the inter-valley 
coupling term is zero (Le., A 2 =0), there are two dou¬ 
ble degenerate energy levels at F point [See left panel of 
Fig. 37 ]. The states at the upper (lower) energy level are 
contributed from “A (B)” sublattices. When ui decreases 
from zero, the inter-valley coupling IA 2 I increases from 
zero, while the sublattice potentials decrease as shown in 
Fig. 3ni As a result, the upper energy level split with a 
gap of 2 A 2 due to the inter-valley scattering, whereas the 
lower one energy level is unchanged since the inter-valley 
scattering only occurs at sublattices “A” as shown in the 
right panel of Fig. 33 As the ui further decrease to 
the critical point, an intersection between |2Ai| and IA 2 I 
indicates that the lower energy level splitting from sub¬ 
lattice “A” reaches the two-fold degenerate energy levels 
of sublattice “B”, for which the bulk band gap closes. 
When ui becomes even smaller, the inter-valley coupling 
is always dominant and the band gap cannot reopen again 
and a chiral single-valley metallic phase is formed. 

METHOD FOR CALCULATING BAND 
STRUCTURE OF PHOTONIC CRYSTAL 

For the simulation of band structure in photonic crys¬ 
tals, we consider a honeycomb lattice with x a/ 3 peri¬ 
odicity, which is comprised of silicon columns linked with 
thin silicon slabs in the vacuum background as shown in 


FIG. 8 : (color online) Left panel: Schematic of honeycomb- 
structured photonic crystals. The distance between nearest 
columns is set to be a and the slab width is d = 0 . 1 a. ri 
(i — 1 - 3 ) label the radii of different columns. Right panel: 
The photonic band structure with uniform column radii and 
the slab widths where Dirac cones are present around K point 
and K' point and the later one is not shown for clarity of the 
figure. 

the left panel of Fig. 3H1 The columns are infinite in 2 :- 
direction and their radii are {i = 1-3). In our calcula¬ 
tion, the electromagnetic wave propagates within the xy 
plane, z.e., the wavevector component along ^-direction 
is kz = 0. Various numerical methods can be used to cal¬ 
culate the photonic band structure, such as plane wave 
method (PWE), finite difference time domain (FDTD) 
method, and finite elements method (FEM) [5|, [6|. Here, 
we used FEM, because it is much efficient in calculat¬ 
ing structures with extremely small domains needing to 
be meshed. The photonic band structure with the uni¬ 
form column radii and the slabs widths are calculated 
as shown in the right panel of Fig. 3H1 where two Dirac 
cones in K and K’ points for transverse-magnetic modes 
[ 3 , 0 are formed, resembling the linear-dispersed Dirac 
cones of pristine graphene. 

HOLLOW AND BRIDGE ADSORPTION 

In above, we study the top-site adsorption case where 
the inversion symmetry is broken and staggered AB sub¬ 
lattice potential is induced. In the following, we show the 
inter-valley coupling mechanisms for hollow and bridge 
adsorption cases with inversion symmetry in 3 x 3 hon¬ 
eycomb super cells. We first consider hollow adsorption 
as shown in Fig. with primitive vectors denoted by 

ai and a 2 in black. By considering only the site ener¬ 
gies surrounding the adatoms, the 7r-orbital in real-space 
tight-binding Hamiltonian is written as: 

iLh = Hq + + ^1^2)5 (10) 

i 

where Yl'i over six atoms nearest to adatoms with 
site energy of ui. Note that AB sublattices are equivalent 
here because of the inversion symmetry. 

The reciprocal lattices for both 1 x 1 (in red) and 
3 X 3 (in black) graphene supercells are displayed in 
Fig. 3H0)- The Brillouin zone of pristine graphene can 
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FIG. 9: (color online) Primitive cells for hollow- (a) and 
bridge-site (c) adsorption, ai indicates the primitive lattice 
vectors for 1 x 1 (in red) and 3 x 3 (in black) graphene su¬ 
percell. ui denotes the site-energy induced by the adatoms, 
(b): Reciprocal lattices for 1 x 1 (in red) and 3 x 3 (in black) 
graphene supercell, (d) and (e): Low energy band structures 
for hollow- and bridge-site adsorptions. A indicates the local 
gap from pseudo-Zeeman held in bridge adsorption. 

be divided into nine copies of that of 3 x 3 graphene su¬ 
percell. By denoting the j-th {j = 1-9) center as Kj, 
the operator can be expanded in momentum space as: 

A' exp[-i(Kj + fe) • where Nq is a 

normalization factor, and k runs over the Brillouin zone 
of 3 X 3 graphene supercell. Therefore, the Hamiltonian 
of Eq. ([TQ|) in momentum space can be expressed as: 

Rh(fe) = Hoik) + ^ (11) 

In the second term, j = j' gives equivalent AB sublattice 
potentials, while j ^ j' couples different parts. Although 
the direct coupling between valleys KK' vanishes due to 
phase interference, Le., ^kk' = = 0, a 

band gap opens at F point with four lower energy bands 
around the gap mainly contributed from eigenstates near 
valleys KK'. As shown in Fig. ^Sjd). This suggests that 
the gap is induced by inter-valley coupling from higher- 
order effects. In below, we demonstrate the physical ori¬ 
gin of inter-valley coupling for the hollow adsorption. 

By doing a block diagonalization @ similar to the top- 
site adsorption case, a low-energy effective Hamiltonian 
at second-order approximation can be reached: 

u 

iJh® = ^ + vpikxCTx + Tzkyay) + -^TxCry, (12) 

o yto 


where r and cr are valley and sublattice Pauli matrices, 
respectively. The first term is an energy shift relative 
to the charge neutrality point, and the second term de¬ 
scribes the kinetic energy with vp being the Fermi ve¬ 
locity. The last term couples valleys K and K', where Tx 
implies a pseudo-Zeeman field in x-direction to induce a 
procession of valley polarization. Moreover, the coupling 
only occurs between different sublattices, and the result¬ 
ing band gap 2u\l^to indicates a second-order correction 
from site energy ui. 

Then we study the bridge adsorption case as shown in 
Fig. 3Hc). Assuming that the adatom only influences 
the site energies ui of the nearest two carbon atoms and 
neglecting the high-order contribution from F valley of 
graphene, the continuum effective Hamiltonian for four 
lower bands can be expressed as follows: 

Hf^ = Y+ Vpikxf^X + TzkyCTy) + (13) 

The first and second terms are the energy shift and ki¬ 
netic energy respectively, whereas the third term rep¬ 
resents the first-order inter-valley coupling contributed 
from the on-site energy which also acts as a pseudo- 
magnetic field in x-direction yet without a sublattice flip¬ 
ping. This term shifts the two degenerate Dirac cones 
of graphene and opens a local energy gap A = 2|rti|/9 
at /c = 0 as shown in Fig. ^Hl^e). The gap is closed 
at {^\ui\/9vf,0) due to the dispersion of energy bands 
where another two Dirac cones are formed. All these re¬ 
sults indicate that the inter-valley coupling mechanisms 
are sensitive to the geometry of adsorption. 


* Correspondence author: qiao@ustc.edu.cn 
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